Draft version November 20, 2012 

Preprint typeset using ETpX style emulateapj v. 5/2/1 1 



ESTIMATING THE "DARK" ENERGY CONTENT OF THE SOLAR CORONA 

Scott W. McIntosh 1,3 , Bart De Pontieu 2 

Draft version November 20, 2012 

ABSTRACT 

The discovery of ubiquitous low-frequency (3-5mHz) Alfvenic waves in the solar chromosphere (with Hin- 
ode/SOT), and corona (with CoMP and the Solar Dynamics Observatory, SDO) has provided some insight into 
the non-thermal energy content of the outer solar atmosphere. However, many questions remain about the true 
magnitude of the energy flux carried by these waves. Here we explore the apparent discrepancy in the resolved 
coronal Alfvenic wave amplitude (~0.5km/s) measured by the Coronal Multi-channel Polarimeter (CoMP) 
compared to those of the Hinode and the Solar Dynamics Observatory (SDO) near the limb (~20km/s). We 
use a blend of observational data and a simple forward model of Alfvenic wave propagation to resolve this 
discrepancy and determine the Alfvenic wave energy content of the corona. Our results indicate that enormous 
line-of-sight superposition within the coarse spatio-temporal sampling of CoMP hides the strong wave flux 
observed by Hinode and SDO and leads to the large non-thermal line broadening observed. While this scenario 
has been assumed in the past, our observations with CoMP of a strong correlation between the non-thermal 
line broadening with the low amplitude, low frequency Alfvenic waves observed in the corona provide the first 
direct evidence of a wave-related non- thermal line broadening. By reconciling the diverse measurements of 
Alfvenic waves we establish large coronal non-thermal linewidths as direct signatures of the hidden, or "dark", 
energy content in the corona, and provide preliminary constraints on the energy content of the wave motions 
observed. 

Subject headings: Sun: chromosphere — Sun: corona — waves 



1. INTRODUCTION 



Alfven waves have been invoked as a driving force behind 
solar coronal heati ng and wind acceleration since their incep- 
tion ([Alfven 1947 ). While Alfvenic waves have been detected 
in situ by several spacecraft in and out of the ecliptic ( Belcher 
1971| ), in the solar atmosphere their presence had only been 



inferred through the la rge non-thermal broadening of coro- 
nal emission lines (e.g.,| Hassler et al. 1990; Hassler & Moran 
T994l|Banerjee et al.||1 998[|Chae, Schuhle & Lemaire| 1998; 
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Banerjee 
Peter I (l999a|b[ |2000a| 
2003; Akiyama, Dosc hek & Mariska||2005] |Singh, Sakurai 
& lcmmoto||2006[ |Singh et al.||2006| ). Only recently were 



they direc tly observe d and measured in the chromosphere ( [De 
Pont ieu et al.|2007a| with the Solar Optical Tel escope (S 



Tsun eta et al.|2008] » on the Hinode spacecraft ( Kosugi et al. 
2007 ) and in the corona with the Coronal Multi-channel Po- 



larimeter (CoM P; ITomczyk et aLl |2007, 2008] |Tomczyk &| 
Mclntosh||2009|). A subsequent investigation with the So- 



lar D ynamics Observator y Atmospheric Imaging Assembly 
(AIA; |Lemen et al.||2011|) confirmed the ubiquity of coronal 
wave motions (Mcintos h et al.|201 1| ). 

Unfortunately, there is a significant discrepancy in the en- 
ergy estimates between these three direct observations. The 
Alfvenic wave energy flux in Doppler motions of the 10747 A 
Fe XIII emission line observed by CoMP is ~0.01 Wm -2 
(Tomczyk et al. 2007 ), some four orders of magnitude less 
than the 100 Win 2 required to balance th e radiative losses of 
the quiet corona ( Withbroe & N oyes|1977] ) and/or accelerate 
the fast solar wind ( Hanstee n & Leer|1995| ). However, the ob- 
servations of SOT indicate that an Alfvenic wave energy flux 
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of order4 kWm~ 2 is present in the chromosphere (De Pontieu 
|et al.|2007a| ) and that even a very low transmission (~3%) of 
these wave s into the corona (due to reflection off the transi- 
tion region, Cranm er & van Ballegooijen 2005 ) will result in 
an energy flux of the order of ^100 Wm z into the corona 
or solar wind. In addition, the presence of at least that much 
wave energy is readily observed in t he quiescent (quiet an d 
coronal hole) corona with SDO/ AIA (Mc intosh et al.|201 1) . 

This leaves us in a quandary: where is the missing wave 
energy in the CoMP observations? Interestingly, the coro- 
nal wave energy flux inferred from the SOT observations (and 
measured by AIA) is consistent with the amplitude o f non- 
thermal broade ning observed by CoMP, as reported by [Tom- 1 
|czyk et alT] ( |2QQ7] ). So, perhaps the real energy flux of Alfvenic 
waves in the corona is hidden from CoMP. In this scenario, 
the vastly different spatial scales of the instruments involved 
would result in considerable superposition of low frequency 
motions resolved by SOT and AIA but unresolved by CoMP, 
quenching the resolved velocity amplitude, energy flux, and 
increasing the non-thermal line broadening of the observed 
line profiles in the optically thin corona. An alternative sce- 
nario faces the challenging task of simultaneously explain- 
ing the large non-thermal broadening in the corona, and the 
observational evidence from AIA of strong Alfvenic waves 
throughout the corona. Perhaps the filling factor of the waves 
observed with AIA is much lower than estimated and/or the 
bulk of the Alfvenic waves (coming up from the chromo- 
sphere) are actually damped and/or dissipated in the first few 
thousand kilometers above the limb, leaving only low re- 
solved wave amplitudes for CoMP? Perhaps the non-thermal 
broadening of the coronal emission line is then caused by un- 
resolved flows into and out of the plane-of-the-sky (POS)? 

In this Letter we investigate these possible scenarios, ad- 
dress the apparent mis-match of Alfvenic wave energies ob- 
served by SOT, AIA, and CoMP and thus constrain the po- 
tential sources of non-thermal broadening, the invariant hid- 
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Fig. 1 . — Comparing the non-thermal linewidth (cr n t) and root-mean- square of the Doppler velocity (Vrms) of CoMP Fe XIII 10774A line with the AIA 193A 
from December 26, 2011. Supporting movies of the CoMP peak intensity, Doppler velocity, and non-thermal width can be found in the electronic edition of the 
Journal. 



den, or "dark" energy, in th e quiet corona. We expand upon 



a sim ple forward model (Mcintosh, De Pontieu & Tarbell 
2008) that explained the center-to-limb variation of transi- 
tion region lines by taking into account that SOT and AIA 
observed a lower atmosphere dominated by a myriad of ex- 
tremely fine-scale, field-aligned and independently moving 
structures, called spicules, with spatial widths <300 km (con- 
siderably smaller than the ^4,500 km CoMP s patial sam- 
pling) that exhibit significant transverse motions ( De Pontieu 
|et aL|2007at|McIntosh et al.|2011) . 

In the following section we discuss the observations that 
drive the model. In Sect.[3| we discuss how the SOT and AIA 
observations are used to drive the model. In Sect. 0] we look 
at the properties of the model that can reproduce the observed 
CoMP behavior. 



2. OBSERVATIONS 
2.1. CoMP Observations 



The CoMP observations used in this paper were obtained 
on December 26, 201 1 from 20: 16 to 21 :45UT in three wave- 
lengths (10745.2, 10746.5, and 10847.8A) across the 10747A 
Fe XIII line, with an exposure time of 250ms at each posi- 
tion. The image groups were obtained at a cadence of 30 
seconds. From the three wavelength positions we measure 
the peak line intensity, Doppler shift, and linewidth of the 
10747 A line. Fitting a histogram to the measured velocities 
allows us to determine the coronal rotation rate and estimate 
the "rest" wavelength to be 10746. 08 A (assuming that half of 
the field of view is red and half is blue- shifted). For an instru- 
mental width (ai ns t) of 21 km/s and a thermal wid th (a t h) of 
21km/s for Fe XIII (using cr t h = y/^ksT* /rrii on , for an ion 
of mass rriion and peak formation temperature of T* assum- 
ing that ion and electron temperat ures are equal) we de termine 
the non-thermal width a nt (= \/ cr l/ e — ^inst ~ G \h)-> wnere 
ai/ e is the 1/e width of the fittea line profile. 

The top row of Fig. [I] compares an AIA 193 A image (panel 
A) taken at the start of the timeseries with the CoMP a nt 
(panel B) and root-mean-square of the Doppler velocity dur- 
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ing the timeseries (Vrms; panel C). The latter is dominated 
by the resolved, low-frequency, and low-amplitude Alfvenic 
waves that are seen to propagate in the corona. Panel D shows 
a two-dimensional histogram of Vrms and cr nt . We find a 
strong correlation between Vrms and a nt . This is an indica- 
tion that the two measures are intimately tied, i.e., the varia- 
tion of the resolved wave amplitude during the timeseries and 
the non-thermal broadening of the line in a pixel/voxel are 
co-dependent. This suggests that whatever causes the strong 
non-thermal line broadening is related to the Alfvenic waves 
that CoMP observes. Note the clear cut-off at cr nt =17-18km/s 
and the enclosing "correlation wedge" of the two dot-dashed 
lines. This wedge is a robust feature of the CoMP observa- 
tions, shifting slightly left and right and tilting slightly up and 
down from day to day as the physical structure of the corona 
above the limb changes. Supporting movies of SDO/AIA 
(Fe XII) 193A, CoMP I, V, and a 1/e can be found in the elec- 
tronic edition of the Journal, the latter clearly illustrate the 
ubiquity of the transverse Doppler motions. 

2.2. Hinode/SOT Observations 

Hinode/SOT observations have shown that the chromo- 
spheric limb is dominated by superposition of a large num- 
ber of long, thin (<200km) spicules that extend to heights of 
5-15" and that undergo significant motions along and perpen- 
dicular to the line-of-sight (LOS). |De Pontieu et alT]pQQ7b| ) 
demonstrated that there are two types of spicules present in 
the Hinode limb observations. Type I spicules evolve on 
timescales of 3-5 minu tes reaching velocities of 10-40km/s 
along their long axis ( |De Pontieu et al.||2007a| ). Type II 
spicules on the other hand, occur on timescales of 10- 120s 
and reach appa r ent long itudinal velocities of 50-150km/s. |De| 
Ponti eu et al.[ (2007b) have shown, using a statistical ap- 
proach, that these chromospheric features undergo vigorous 
Alfvenic motions with amplitudes around 20(±5)km/s and 
periods ranging from 100 to 500s. 

2.3. SDO/AIA Observations 



Mcintos h et al. ( 2011| ), using a similar approach, demon- 
strated that an energetically relevant portion of the Alfvenic 
energy observed in the chromosphere by SOT is indeed visi- 
ble in emission formed at coronal temperatures. Those waves 
have made it "into" the corona - these observations were made 
possible by the high S/N, low scattered light, high spatial res- 
olution observations of AIA. The imaged wave motions have 
amplitudes, periods, and speeds commensurate with those 
observed by SOT near the limb, and can b e considered the 
coronal extension ( |De Pontieu et al."||2011|) of the ch romo- 
spheric features studied by |De Pontieu et al.] ( |2007b| ). The 
observed phase speeds of these coronal Alfvenic waves reach 
500-1000km/s, values that are commensurate with those ob- 



served by CoMP dTomczyk et al. 20071 ITomczyk & Mcintosh 
|2U09] ). 

3. METHOD & ANALYSIS 

We use the high resolution SO T ( |De Pontieu et al.|[ 2007a) 
and AIA (Mcintosh et al. 2 01 1) observations as a guide for 
the transverse motions that occur on a spatial scale that is far 
smaller than CoMP resolution. We develop a (simple) Monte- 
Carlo forward model of Alfvenic waves in the corona based 
on the concept of a "thread", or elementary oscillating struc- 
ture. Our simple picture is that any pixel in the optically thin 
model corona can be populated with N independent threads 



and the 10747 A signal observed is given (at any time) by sum- 
ming the emission of all of the threads in any given pixel. 

The properties of any single thread are as follows: a thread 
can appear at any time with a lifetime chosen from a Gaus- 
sian distribution of 100 ± 20s; has a uniform brightness; and 
is subject to Alfvenic motion governed by a randomly chosen 
period, amplitude and phase. Here we consider swaying mo- 
tions, since those are clearly observed with AIA. We assume 
a polarization angle chosen from a uniform distribution be- 
tween and 2tt with respect to the LOS, and apply a 1/ cos 6 
correction factor to accommodate the fact that not all threads 
are radial (where follows a Gaussian distribution of 0±20°). 
The sinusoidal Alfvenic motions of the thread have a period 
drawn from a distributi on based on the CoM P power spec- 
trum (see, e.g., Fig. 2 of |Tomczyk et al.|2007| ), a phase from a 
uniform distribution between and 2n and an amplitude from 
a Gaussian distribution around X±5 km/s (where X can be 
any value >5km/s). Each thread emits a Gaussian line profile 
of width 21 km/s (equivalent to a t h of the Fe XIII plasma) at 
the given LOS velocity determined from its period, amplitude 
and phase. The resulting line profile in a pixel is the sum of 
the Gaussians of the N individual threads. 

We then observe these sy nthetic profiles in the same fashion 
as the CoMP observations ( [Tomczyk et al.|200"7] ) by sampling 
the 10747 A line at the three wavelengths above with a Gaus- 
sian filter width of 1.3 A, and fitting the ensemble line profiles 
with a single Gaussian to derive I, V, and cri/ e . Figure|2]shows 
sample (pixel) timeseries of the Doppler velocity (ancTV rms ; 
left column) and non-thermal linewidth (<r nt ; right column) 
for a range of thread numbers (bottom to top: 200, 500, 1000, 
and 2000) and a fixed 40km/s wave (peak-to-peak) ampli- 
tude (Vp2p). Notice the diminishing values of Vrms with 
increasing thread number - these values eventually saturate 
at the Gaussian fitting accuracy of the three-point measure- 
ments (~0.1 of the spectral resolution). Notice also that a nt 
& >/2Vp2p and the gradual increase in <j nt with increasing 
thread number and the reduction of its variance with time - 
recalling that Tom czyk et al.| ( 2007 ) noticed an invariance of 
the linewidths with time in their original measurements. To 
investigate any artifacts introduced by the three-point mea- 
surements (red in Fig. [2]), we also perform the same calcu- 
lations of velocity and linewidth using synthetic data at con- 
siderably higher spectral resolution (blue). We find that there 
is no significant difference in Vrms or cr nt for both spectral 
resolutions except in the case of large thread numbers where 
the profiles can become non-Gaussian. In summary, it is clear 
from Fig. [2] that our model can reproduce both the low am- 
plitudes ofthe resolved Alfvenic waves and large temporally 
invariant non-thermal line broadening observed with CoMP. 

Extending this computation to a broader range of models 
is straightforward. The panels of Fig. [3] show the variation in 
Vrms (A) and ai/ e (B) for a range of thread numbers and 
input wave amplitudes where the other thread parameters re- 
main as described above. The white (dot-dashed and solid) 
lines show the 15 and 45km/s values of cr nt respectively, while 
the black (dot-dashed and solid) lines correspond to the 0.25 
and 1.5km/s values of Vrms respectively. These bounds en- 
capsulate more than 75% of the analyzed pixels in panels B 
and C of Fig. [T] Note that iso-contours of Vrms are approx- 
imately quadratic functions of the (log of the) thread number 
and input wave amplitude, whereas the measured linewidths 
almost directly reflect the input wave amplitude of the model 
(for thread numbers larger than ~30). 

One of the most stringent constraints from CoMP for our 
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Fig. 2. — Sample timeseries for different numbers of emitting threads in a coronal voxel (increasing from bottom to top: 200, 500, 1000, and 2000) for the 
three-point CoMP measurements (red) and an imaginary spectrograph with 50 times the spectral sampling (blue). The left column of panels show the variation 
of the Doppler velocity and Vrms (green dot-dashed line) while the right column shows the variation of cr n t for a sample timeseries of 200 minute duration. 



model is the correlation and peculiarly shaped correlation 
wedge that we found in Fig.[Tp between the Vrms and <r nt . 
We show in Fig. H] that we can easily (and naturally) repro- 
duce such a correlation of Vrms and a nt and the wedge by 
using the data from the simple model presented in Fig. [5] We 
use the line-enclosed regions in the panels of Fig [3] (which 
are driven by CoMP measurements), and apply other simple 
selection rules on the parameter set to best recover the corre- 
lation wedge. We find that the shape of the wedge is sensitive 
to the number of threads along the line of sight, the mini- 
mum wave amplitude and any additional broadening not due 
to swaying motions. Panel A shows the relationship between 
Vrms and a nt when the minimum wave amplitude is allowed 



to be zero. To match the cut-off of the wedge at 17km/s shown 
in the observations of Fig. [T] we need to add 17km/s of addi- 
tional broadening. The maximum and minimum number of 
threads determines the lower and upper slopes of the wedge. 
While the slopes of the wedge are generally well reproduced, 
the low end of the wedge is too pointed and the scatter in 
Vrms does not fill the wedge. This implies that this param- 
eter set is inadequate to reproduce the observed relationship. 
Similarly, in panel B, we have omitted the additional broad- 
ening but, to match the observed low-end cut-off at 17km/s 
we must impose a minimum wave amplitude of 23km/s peak- 
to-peak (which is equivalent to 16km/s rms). This parameter 
set is also inadequate: the slope of the wedge is too shallow 
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Fig. 3 . — Variation of Vrms (A) and a 1 / e (B) measured from the Monte- 
Carlo simulated timeseries of three-point 10747 A measurements for a range 
of thread numbers and input wave amplitudes. The white (dot-dashed and 
solid) lines show the 15 and 45km/s values of cr n t respectively, while the 
black (dot-dashed and solid) lines correspond to the 0.5 and 1.5km/s values 
of Vrms respectively. 



and it does not fill the wedge at values of high Vrms and 
cr nt . Finally, in panel C, we show the optimal parameter set 
for this observation. The best match fills the wedge and re- 
quires a minimum wave amplitude of 14km/s peak-to-peak 
(which is equivalent to lOkm/s rms) and an additional broad- 
ening of 14km/s rms. The high a nt cut-off would indicate that 
the maximum wave amplitude present is 40km/s (or 56km/s 
peak-to-peak). The upper bound of the wedge is determined 
by pixels that have low thread numbers where we find the 
best match is at 125 threads per voxel. The lower bound of 
the wedge is determined by the highest number of threads in a 
voxel where the best match is between 1200 and 1300 threads. 
In summary, the correlation found from CoMP comes about 
naturally if the observations are dominated by superposition 
of vigorous Alfvenic motions (with properties similar to SOT 
and AIA) along a LOS with a large and variable number of 
independently oscillating structures. 

4. DISCUSSION 

To understand the low values of the resolved CoMP wave 
amplitude we have constructed a simple Monte-Carlo model 
with only two variables: the input wave amplitude, and the 
number of emitting threads in a voxel. Our model can repro- 
duce the CoMP observations, even though it is based on the 
much larger wave amplitudes observed with SOT and AIA: 
CoMP measurements will see decreasing wave amplitudes 
and increasing linewidths with increasing complexity (super- 
position) of the emission along the LOS. We also see that the 
linewidths approach temporal invariance within the accuracy 
of the Gau ssian fit - this invariance was noted by (Tomczyk 
et al.|2007j ). We note that the movies supporting Fig. [T] show 
that in the less dense regions (in coronal holes), the linewidths 
do show some periodicity with time. This would appear to be 
consistent with the lower panels of Fig. [2] which show low am- 
plitude variations of linewidths where the number of threads 
is sm aller. The earlier dedu ction of complete temporal invari- 
ance ( [Tomczyk et al. 2007 ) may have been biased by the study 
of a single sub-region. 

We have also reproduced the wedge relationship between 
Vrms and a nt . This relationship comes about naturally if 
one limits the input wave amplitudes and assumes a min- 
imum and maximum number of independently oscillating 
threads along the LOS. Under the assumptions of our sim- 
plistic model, this would imply that observed CoMP wave 
amplitudes of ~lkm/s are actually hiding the true wave am- 
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Fig. 4. — Scatter plots of the Vrms and cr nt values used to reproduce 
the variance observed in panel D of Fig. [T] The red dashed lines are taken 
from the bounds of the correlation wedge in Fig.ffl Panel A has no minimum 
wave amplitude, minimum and maximum threacT numbers of 110 and 1200 
(respectively), but requires 17km/s of additional broadening to match the ob- 
served cut-off velocity of the wedge. Panel B has non additional broadening, 
minimum and maximum thread numbers of 300 and 1500 (respectively), but 
requires 23km/s minimum wave amplitude to match the observed cut-off ve- 
locity of the wedge. Panel C is the optimal solution with 14km/s additional 
broadening, a 14km/s minimum wave amplitude and thread limits as shown. 



plitudes that range between 25 and 56km/s. Further, the mea- 
sured linewidths show an almost one-to-one relationship with 
the amplitude of the waves present (FigJSb), thus providing, 
in principle, an almost direct measure ofthe unresolved low- 
frequency wave energy. We note that the exact calibration 
between line broadening and wave amplitude of swaying mo- 
tions is actually dependent on several assumptions (see be- 
low). 

Unresolved flows along the LOS likely contribute to the 
non-thermal line broadening observed with CoMP. In prin- 
ciple, such flows could even cause the correlation found in 
Fig.[Tp, as long as the flows are temporally varying (Vrms)- 
However, this is unlikely for a couple of reas ons. First, the 
Vrms is not dominated by flows but by waves ( [Tomczyk et al.| 
2007). Second, as shown in Fig. [5] the upward-downward 
directed wave power ratio of (Tomczyk & Mcintosh 2009) 
shows the same correlation with the cr nt as Vrms- Dividing 
the coronal cavity into its southern (foreground) and north- 
ern (background) portions by the dashed radial line we plot 
the correlation of a nt and the directional wave power ratio. 
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In both regions the relationship is linear, which supports the 
scenario in which cr nt is dominated by the propagating wave 
motions of the plasma and not by bulk coronal flows along the 
LOS. 

Our simplistic model has several limitations. For exam- 
ple, we assume that the transverse motions are dominated by 
swaying, and not torsional, motions. The latter would not 
cause resolved Doppler shift variations with time. It is pos- 
sible that the additional broadening of 14km/s (uncorrelated 
with Vrms) required to reproduce the correlation wedge is 
a signature of such torsional motions. This would be con- 
sistent with recent observations of peak-to-peak ampli tudes 
in spicules of order 20-30 km/s ( |De Pontieu et al.|[20T2b . It 
is also possible, as recent studies suggest, that some of the 



Fig. 5. — Compari ng cr n t (panel A ) and directi onal wave power ratio (panel B) from CoMP in the coronal cavity where the wave signal to noise is the highest 
as studied by Tomczyk et al. ( 2007) and Tomczyk & Mcintosh (2009). The cavity is divided into a North and South component by the dashed radial line. The 
right panel shows the correspondence between these quantities for both portions of the coronal cavity. 

propagate upward from the loop footpoints. 

Taken at face value, our model would predict a very large 
number of independently oscillating threads along the LOS 
(~1 solar radius) contained in one CoMP pixel (4,500 km). 
It may well be the case that the corona is extremely finely 
structured. However, the number of threads predicted by our 
model may be very different from the number of threads in 
the corona. This is because the exact number will critically 
depend on the mix of torsional and swaying motions, and on 
the (perpendicular) coherence of the waves. The latter is un- 
known: if phase mixing plays a significant role, wave am- 
plitudes seen by CoMP can, in principle, be quenched much 
more efficiently than the completely independently oscillating 
threads assumed here. 

Our data seem to exclude a scenario whereby high fre- 
quency Alfvenic waves (invoked by several models for the so- 
lar wind) dominate non-thermal line broadening (at least for 
heights <50Mm). If the wave generation process low down 
predominantly created high frequency waves, with a small 
amount of wave power at low frequencies, one could in prin- 
ciple recover the correlation wedge. However, the AIA and 
SOT observations of strong, low frequency Alfvenic swaying 
at low heights argue against such a scenario. 

NCAR is sponsored by the National Science Founda- 
tion^ CoMP data can be found at |http : //mlso . hao . | 
ucar.edu/ Hinode is a Japanese mission developed and 
launched by ISAS/JAXA, with NAOJ as a domestic partner 
and NASA and STFC (UK) as international partners. We 
acknowledge support from NASA contracts NNX08BA99G, 
NNX11AN98G, NNM12AB40P, NNG09FA40C (IRIS), and 
NNM07AA01C (Hinode). 



swaying motions couple to torsional modes ([Verth, Terradas, 
& Goossens|2010l|Pascoe et al.|2011HDe Moortel & Pascoe 



2012| ). If such a co upling mechanis m occurs in th e corona 



(see, e.g., |Ionson|1 978; Rud erman & Roberts |2002| ), the am- 
plitude of the torsional modes could be correlated with that 
of the swaying motions, so that the line broadening caused by 
torsional modes could actually be part of the broadening that 
is correlated with Vrms- In this case, torsional motions would 
not cause the additional broadening (which instead could be 
caused by unresolved flows along the LOS). If this were the 
case, the observed correlation between Vrms (from swaying) 
and a i / e would not necessarily all be caused by LOS superpo- 
sition. In this scenario, the relative contributions of swaying 
and torsional motions likely change with height (with increas- 
ing contributions from torsional motio ns higher up), because 



the co upling mechanism suggested by |De Moortel & Pascoe 
( 2012) occurs over a finite time /distance during which waves 



Alfven, H., 1947, MNRAS, 107, 211 
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